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Abstract: Osteoporotic fracture carries an enormous public health burden in terms of mortality 
and morbidity. Current approaches to identify individuals at high risk for fracture are based on 
assessment of bone mineral density and presence of other osteoporosis risk factors. Bone mineral 
density and susceptibility to osteoporotic fractures are highly heritable, and over 60 loci have been 
robustly associated with one or both traits through genome-wide association studies carried out 
over the past 7 years. In this review, we discuss opportunities and challenges for incorporating 
these genetic discoveries into strategies to prevent osteoporotic fracture and translating new 
insights obtained from these discoveries into development of new therapeutic targets. 
Keywords: bone mineral density, genome-wide association studies, osteoporosis, prediction, 
fracture, genetics 

Introduction 

Osteoporosis is a major public health problem that is associated with significant mortal- 
ity and morbidity. Although osteoporosis is often regarded as a disorder of women, and 
indeed women are disproportionately affected the disorder also carries an enormous 
health burden in men. An estimated 20% of women and 4% of men in the US were 
diagnosed with osteoporosis in 2005-2006, with rates increasing to 48% in women 
and 12% in men over the age of 85 years. 1 It has further been estimated that up to one 
half of all women and one quarter to one third of all men over the age of 50 years will 
experience a hip fracture during their lifetime. 2 The economic impact is enormous, 
with an estimated health care cost of $17 billion in direct costs in the US alone. 3 

Osteoporosis is a systemic skeletal disorder characterized by microarchitectural 
deterioration in bone tissue that results in reduced bone strength, bone fragility, and 
increased susceptibility to fracture. 4 The most common fracture sites are the hip and 
spine, although any bone can be affected. Identifying individuals with low bone 
strength and high risk for fracture is important so that strategies to slow bone loss, 
including pharmacotherapy, can be deployed. Although the molecular basis underly- 
ing these defects and the triggers that initiate them are not completely understood 
a strong familial contribution to bone mineral density (BMD) is well established with 
genes estimated to contribute 60%-80% of the variability in BMD. 5 6 It is hoped that 
current efforts aimed at identifying the specific genes involved will provide further 
insight into the molecular basis underlying osteoporosis and will eventually lead to 
new therapeutic targets and prevention strategies. 

The goal of this review is to summarize the current knowledge of factors contribut- 
ing to increased fracture risk in osteoporosis, with an emphasis on genetic factors, and 
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to discuss the potential impact of new genetic discoveries on 
identifying individuals at risk for osteoporotic fracture and 
preventing future fractures. 

Development and sequelae 
of osteoporosis 

Bone is a dynamic tissue that is continually being formed and 
resorbed during one's lifetime. Bone mass at any time of life 
is the integral of the amount of bone accrued during growth 
and consolidation (peak bone mass) and the loss of bone that 
occurs with aging. The importance of achieving an optimal 
level of bone mass in early adulthood has highlighted the role 
of early nutritional and lifestyle factors in bone health, and 
led some to characterize osteoporosis as a pediatric disease, 7 
or at least an adult disease with its seeds in childhood. 

There is substantial variability both in the acquisition of 
peak bone mass, attained during early adulthood (by women 
in the early 20s and by men in the mid 20s), and in the rate 
of bone loss that begins during the perimenopausal years. 
In both men and women, peak bone mass is generally main- 
tained in the 30s and early 40s, with a slow loss of bone mass 
that accelerates in women during the menopause. Both peak 
mass and rate of bone loss contribute to fracture risk; that is, 
the fracture risk will be high in individuals who achieved a 
low peak bone mass even if the subsequent rate of bone loss 
is small, as well as in individuals who achieve a relatively 
high peak bone mass but who experience a relatively high 
rate of bone loss in their later years. In a given year, about 
10% of the skeleton is broken down and reformed by bone 
remodeling. 8 

Bone remodeling occurs at numerous focal areas throughout 
the skeleton. The net amount of bone formation and resorp- 
tion can be crudely assessed by bone turnover markers that 
are routinely measured in many laboratories. Bone formation 
markers, which are protein products of the osteoblast, include 
bone-specific alkaline phosphatase, procollagen- 1 amino 
peptide, and osteocalcin, while resorption markers include 
C-telopeptide and N-telopeptide, which are breakdown prod- 
ucts of collagen (the main protein in bone). 

Osteoporosis is associated with increased morbidity and 
mortality. In the first 6 months following hip fracture, there is 
an excess mortality of at least 18%.' A meta-analysis showed 
that mortality 1 year after hip fracture was more than twice as 
high in men as in women. 10 Although men are younger than 
women at the time of hip fracture, they fare worse than women 
following hip fracture due to greater comorbidity and higher 
postoperative infection rates (reviewed by Sterling 11 ). Increased 
mortality also follows vertebral compression fracture. 12 



In addition to increased mortality, osteoporotic fracture is 
frequently followed by chronic pain, functional limitation, 
and psychologic difficulties, that include fear of falling and 
social isolation. 13 Overall quality of life is reduced in both 
women 14 and men 15 following osteoporotic fracture. Medical 
treatment of osteoporosis reduces fracture risk and enhances 
quality of life, 16 and some medications have been associated 
with improved overall survival. 17 

Diagnosis, risk estimation, 
and treatment 

Bone strength is determined by two factors, ie, BMD and 
bone quality. 18 BMD is determined by dual energy X-ray 
absorptiometry (DXA) and is usually performed at the 
lumbar spine and hip. Bone quality can be assessed by three- 
dimensional imaging modalities, such as peripheral computed 
tomography and high resolution magnetic resonance imag- 
ing, 19 although these are currently used primarily in research 
settings and are not available clinically. The assessment of 
fracture risk using only BMD has limitations. 

In postmenopausal women and men over the age of 
50 years, the diagnosis of osteoporosis can be made either 
by the occurrence of a fragility fracture (following a fall 
from standing height, after ruling out malignancy) or by 
a DXA T-score of ^—2.5 at the lumbar spine, total hip, or 
femoral neck (Table 1). This T-score threshold corresponds 
to a BMD that is 2.5 standard deviation (SD) units or more 
lower than the BMD for an "average" young adult female. 
The World Health Organization has published guidelines 
further classifying a T-score £—1.0 as normal and a T-score 
of -1.0 to -2.5 in the absence of fracture as "low bone 
mass" (or osteopenia). 20 The definition of osteoporosis in 
premenopausal women and men younger than 50 years is 



Table 1 Diagnosis of osteoporosis 




Postmenopausal W 


Premenopausal W 




and M >50 years 


and M <50 years 


Osteoporosis 


T-score £—2.5 or 


Z-score £—2.0 and 




fragility fracture 


fragility fracture 


Low bone mass 


-1.0 to -2.5, 




(osteopenia) 


no fracture 




Low bone mass 




Z-score <— 2.0 and 


for age 




no fracture 


Normal 


T-score £—1.0, 


Z-score s— 2.0, no 




no fracture 


fracture 



Note: Adapted from: World Health Organization. Assessment of fracture risk and 
its application to screening for postmenopausal osteoporosis. Report of a WHO 
Study Group. World Health Organ Tech Rep Ser. 1994;843:1-129 20 and Martinez- 
Morillo M, Grados D, Holgado S. Premenopausal osteoporosis: how to treat? 
Reumatol Clin. 20 1 2;8:93-97. 21 
Abbreviations: W, women; M, men. 
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controversial, but most agree that the patient must have expe- 
rienced a fragility fracture and have a BMD <2 SD that of 
the average BMD for individuals of the same age and gender 
(ie, a Z-score ^-2.0). 21 

Pharmacologic treatment to prevent fracture is cur- 
rently targeted to those at high risk of fracture. According 
to guidelines established by the US National Osteoporosis 
Foundation, postmenopausal women and men older than 
50 years and with osteoporosis (fracture orT-score ^-2.5 at 
the spine or femoral neck) should be treated. 22 Current phar- 
macotherapies to prevent osteoporotic fracture are designed 
to slow bone loss. A number of different drugs are available, 
most of which act as antiresorptive agents on osteoclasts to 
decrease bone loss, although there is currently one that acts 
by stimulating bone formation (teriparatide). There is a wide 
variation among patients in response to drugs, both in terms 
of efficacy and frequency/type of adverse reactions. 23 ' 24 

Because of the high prevalence of osteoporosis and 
the proven benefit of pharmacologic treatment in reducing 
fractures, 25 there is a strong incentive to identify individuals at 
risk for future fracture. DXA is the most widely used tool for 
this purpose, although there are concerns that DXA utilization 
rates have dropped since 2007, coincident with a lowering 
of Medicare reimbursement rates for non-facility-based 
DXAs. 26 " 28 Nevertheless, in a large meta-analysis involving 
approximately 90,000 years of follow-up and over 2,000 
fractures, it was estimated that each SD decrease in BMD 
from the age-adjusted mean was associated with a 2.3-fold 
increase in incidence of vertebral fractures and a 2.6-fold 
increase in hip fractures. The predictive ability was roughly 
similar to (or, for hip or spine measurements, better than) that 
of a one SD increase in blood pressure for stroke and better 
than a one SD increase in serum cholesterol concentration 
for cardiovascular disease. 29 In a later meta-analysis, Johnell 
et al reported that each SD decrease in BMD corresponded 
approximately to a 2.9-fold increase in hip fracture risk. 30 

While these results illustrate the utility of DXA as a 
screening tool, DXA does not have a high sensitivity at the 
individual level for fracture prediction. 31 - 32 Indeed while the 
lifetime risk of fracture in women and men over age 50 years 
is estimated to be 50% and 30%, respectively, 2 more than 50% 
of women and 70% of men who experienced a fracture did not 
meet DXA criteria for osteoporosis. 33 Tools that assess bone 
quality (eg, high resolution peripheral computed tomography) 
improve fracture risk prediction over DXA but are not avail- 
able clinically. 34 To improve fracture prediction, the FRAX® 
tool ( http://www.shef.ac.uk/FPvAX/ ) was developed in 2008 
by the World Health Organization as an online calculator that 



integrates information from BMD and clinical risk factors 
to calculate the 10-year probability of hip fracture and the 
10-year probability of any major osteoporotic fracture (clini- 
cal spine, forearm, hip, or shoulder fracture). The clinical 
risk factors included in FRAX fracture prediction include 
age, gender, height, weight, prior fracture, parental history 
of hip fracture, current glucocorticoid use, current smoking, 
excess alcohol intake, rheumatoid arthritis, and secondary 
causes of osteoporosis. Separate prediction equations have 
been developed for different geographic regions and ethnic 
groups. The National Osteoporosis Foundation guidelines 
recommend pharmacologic treatment of those with low bone 
mass who have a 10-year risk of fracture of at least 20% for 
any fracture or 3% for hip fracture as calculated by FRAX. 22 
Using these treatment guidelines, approximately 30% of 
women and 19% of men over age 50 years in the US meet 
the criteria for pharmacologic treatment. 35 

Risk factors for osteoporotic 
fracture 

Because of the wide availability of DXA, the high precision 
with which it can be measured and the correlation of BMD 
with fracture risk, much of the research into the determinants 
of osteoporotic fracture has focused on identifying the risk 
factors associated with variability in BMD. We provide 
below and in Table 2 a brief summary of risk factors for 
BMD and fracture; more extensive reviews on this subject 
have been published elsewhere. 36-38 

BMD is lower and hip fracture rates are higher in women 
compared with men and in populations of Caucasian origin 
compared with non-Caucasian populations. Higher body 
weight is also consistently associated with higher BMD and 
lower hip fracture risk, a correlation that may be attributable 
to multiple factors, including higher mechanical loading on 
skeletal bone that may promote mineralization and bone 
strength, skeletoprotective effects of adipokine secretion 
from fat tissue, and cushioning that may protect bone from 
the trauma of falls. Most non-vertebral fractures are associ- 
ated with falls, and this dependence on falls complicates the 
ability to predict fractures accurately. 

Reproductive and hormonal factors play an important role 
in regulating bone metabolism, with bone loss accelerating in 
women at the onset of menopause. Bone sites most affected 
by estrogen deficiency are those with high trabecular content, 
such as the spine and ultradistal forearm. 36 Estradiol is the 
most important hormone for bone health in both genders, but 
other hormones are important, including growth hormone, 
testosterone, and insulin. Nutritional factors, including calcium, 



The Application of Clinical Genetics 2013:6 



submit your manuscript | www.dovepress.com 
Dovepress 



77 



Mitchell and Streeten 



Dovepress 



Table 2 Non-genetic risk factors for osteoporosis 



Personal factors 



Medications 



Diseases 



Caucasian 
Female 

Age >S0 years 
Low body weight 

Physical inactivity 
Falls 

Muscle weakness 
Reduced vision 
Smoking 
Excess alcohol 

Elevated homocysteine (females only) 



Glucocorticoids 

Anticoagulants (warfarin, heparin) 

Anti epileptics 

Sex hormone suppressants 

(eg, aromatase inhibitors, Depo-Provera®) 

Proton pump inhibitors 

Thiazolidinediones 

Antidepressants 



Monogenic diseases 

Ol, OPPG 
Muscular dystrophy 

Connective tissue disorders (eg, Marfan syndrome, 
Ehlers-Danlos syndrome) 
Cystic fibrosis 

Metabolic (eg, glycogen storage diseases) 
Complex diseases 

Renal (chronic kidney disease, hypercalciuria) 
Rheumatoid arthritis 
Hematologic (eg, mastocytosis) 
Gastrointestinal (inflammatory bowel disease, 
celiac disease) 

Pulmonary (chronic obstructive pulmonary disease) 
Endocrine (types I and 2 diabetes mellitus, 
hyperthyroidism, hyperparathyroidism) 



Abbreviations: Ol, osteogenesis imperfecta; OPPG, osteogenesis imperfecta ocular form. 



vitamin D, and protein intake, play an important role in both the 
acquisition of peak bone mass and its maintenance in later life. 
However, a more precise understanding of the role of nutritional 
factors on BMD and bone health has been hampered by meth- 
odologic limitations in assessing dietary intakes and the time 
lag between dietary intake and its effects on bone health. 

Other important risk factors for osteoporosis and reduced 
BMD include physical inactivity (which may increase frac- 
ture risk in part by reducing muscle mass) cigarette smoking, 
alcohol use, low sun exposure (for production of vitamin D), 
and use of some medications, including glucocorticoids and 
anticonvulsants. The impact of some of these osteoporosis 
risk factors differs slightly according to age and across sites, 
perhaps due in part to differing compositions of cortical 
versus trabecular bone. 39 Recent studies have also indicated 
that hip geometry contributes to fracture risk. 40,41 

Recent studies of genetic factors 
associated with osteoporosis 

There is a strong genetic contribution to susceptibility to osteo- 
porosis, with genes estimated to account for about 25% of the 
variance in liability to osteoporotic fractures, 42 25%-54% for 
fractures of the wrist, 43 - 44 and up to 48% for fractures of the 
hip. 42 Despite these moderate to high heritabilities, the fracture 
phenotype is a very challenging one for use in genetic stud- 
ies because fracture risk is influenced by a number of diverse 
physiologic factors, including BMD and age-related declines 
in bone microarchitecture and quality, muscle strength, bal- 
ance, cognition, cardiovascular function, and vitamin D status. 
Since each of these factors is itself under at least partial genetic 
control, variants that influence fracture susceptibility entirely 



through any of these other factors should be more easily detect- 
able in an analysis of the factor itself, rather than fracture. 

Because of its wide use and acceptability as a screening 
tool, numerous genetic studies have been carried out on BMD, 
providing highly consistent estimates that genes account for 
60%-80% of the total variability in BMD. 5 That low BMD 
is a significant mechanism through which fracture risk can 
be transmitted through families was demonstrated by See- 
man et al over 20 years ago, who showed that BMD in the 
lumbar spine and femoral neck was lower in the daughters 
of women with osteoporotic fractures than in the daughters 
of non-osteoporotic women. 45 Significant heritabilities of 
BMD in young adulthood and premenopausal adults have 
been reported, 46 ' 48 implying a genetic contribution to the 
acquisition of peak bone mass. It is also presumed that genes 
contribute significantly to variability in aging-related bone 
loss, and in fact a few small studies have reported moderate 
heritabilities of bone loss. 46 49 " 52 However, precise estimates 
of the genetic contribution to bone loss are lacking because 
of the difficulty in assessing bone change reliably, ie, specific 
trajectories of bone change are non-linear and highly age- 
dependent, making it critical that multiple measurements 
over time be obtained in a standardized fashion. 53 These dif- 
ficulties notwithstanding, it would be immensely valuable to 
identify bone change-related genes because many may fall in 
pathways that are more amenable to therapeutic modification 
than pathways specific for bone acquisition. 

Early efforts to identify specific genes related to variation 
in BMD and fracture risk focused on identifying biologically 
motivated candidate genes and testing specific genotyped 
variants for association with BMD (or fracture). While many 
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positive association results were published, with few excep- 
tions (eg, the estrogen receptor 1 [ER1 gene] and low-density 
lipoprotein receptor-related proteins 4 and 5 [LRP4, LRP5] 
genes), most initially reported associations turned out to be 
difficult to replicate. With advances in genomic technology, 
numerous genome-wide association studies (GWAS) of BMD 
and related traits have been published over the last 6-7 years. 
It is important to appreciate that the GWAS approach is 
designed to test the hypothesis of association with genetic 
variants that are common in the population, ie, typically vari- 
ants with minor allele frequencies of 5% or greater. Although 
the variants (typically single nucleotide polymorphisms 
[SNPs]) included on GWAS arrays are generally neutral (ie, 
do not have functional effects), they may be highly correlated 
with numerous other surrounding SNPs, including some 
that may influence gene expression or function. The largest 
of the BMD GWAS studies published to date was from the 
Genetic Factors of Osteoporosis Consortium and included a 
discovery sample of over 32,000 subjects from 17 different 
populations and an independent replication sample of over 
50,000 subjects. 54 In this study, 56 loci significantly associ- 
ated with BMD, of which 14 were found to associate also 
with fracture. Several additional BMD-associated loci have 
since been identified, some from non-Caucasian populations, 
bringing the total number of BMD-associated loci to 62. 55 The 
loci that have been robustly associated with BMD through 
GWAS studies to date are listed in Table 3. The effect sizes 
of these loci are uniformly small, each accounting for far less 
than 1% of the total variation in BMD; even collectively, all 
identified loci account for less than 6% of the total variation 
in femoral neck BMD. 54 Future GWAS studies including even 
more subjects will undoubtedly be published, through which 
even more BMD-associated loci will be identified. However, 
these newly associated loci will almost certainly have even 
smaller effect sizes than those already discovered. 

One goal of GWAS is to obtain insights into disease 
pathogenesis by determining whether the associated loci map 
to any novel molecular pathways. Many of the 62 associated 
loci do indeed map to genes that fall within pathways related 
to bone biology, such as the Wnt/(3-catenin signaling path- 
way, the RANK-RANKL-osteoprotegerin (OPG) pathway, 
mesenchymal stem cell differentiation, and endochondral 
ossification, although the involvement of these pathways 
in osteoporosis pathophysiology had generally been known 
before the advent of GWAS. The Wnt/(3-catenin signaling 
pathway, which is mediated by interactions of Wnt proteins 
with their receptors, causes an accumulation of (3-catenin in 
the cytoplasm that then translocates into the nucleus where 



it participates in gene transcription. The importance of this 
pathway in bone health was first recognized with the discovery 
of a mutation in LRP5 as the cause of osteoporosis pseudo- 
glioma syndrome, 56 and it is now appreciated that this pathway 
is crucially important for a variety of processes, including 
bone cell differentiation, proliferation, and apoptosis. The 
RANK-PsANKL-OPG pathway is an important regulator of 
bone resorption that involves receptor activator of nuclear 
factor-KB (RANK), its ligand (RANKL), and OPG, a so- 
called decoy receptor of RANKL. RANK is expressed by 
osteoclasts and their precursors, RANKL is expressed on 
osteoblast surfaces, and OPG is produced by osteoblasts. It is 
the binding of RANKL to its receptor, RANK, that controls 
the differentiation, proliferation, and survival of osteoclasts. 
Figure 1 illustrates how the Wnt and RANK/RANKL/OPG 
pathways interact with each other to regulate the balance 
between bone formation and resorption. 

The process of mesenchymal stem cell differentiation is 
highly relevant to bone turnover because mesenchymal stem 
cells are multipotent stromal cells that can differentiate into a 
variety of cell types, including osteoblasts, chondrocytes, and 
adipocytes. The endochondral ossification pathway involves 
processes active during fetal development of the mammalian 
skeleton. A fuller description of the Wnt, RANK-RANKL- 
OPG, and endochondral ossification pathways is beyond the 
scope of this review, but excellent reviews are available on 
these subjects. 57 

Finding genes for fracture risk is likely to be more dif- 
ficult than for BMD due to the much smaller sample sizes 
generally available for studies of fracture and the complexity 
of the fracture phenotype. As previously noted, it is difficult 
enough to identify genes/SNPs associated with intermediate 
traits. For example, although over 60 SNPs have now been 
associated with BMD, for which the heritability is very high, 
the effect sizes of all are very small, and enormous sample 
sizes were required to identify these. Of the 16 SNPs that 
have been associated with fracture to date (see Table 3), most 
were tested because of their initial association with BMD, 54,55 
and all have odds ratios for fracture of 1.11 or lower for 
the risk allele with the exception of one, ie, rsl3 182402 in 
ALDH7A1 (odds ratio 2.25). This SNP was identified in a 
GWAS of fractures in a Chinese population. 58 ALDH7A1 is 
a gene in the aldehyde dehydrogenase 7 family (member Al) 
that degrades and detoxifies acetaldehyde, which inhibits 
osteoblast proliferation and results in decreased bone forma- 
tion. 59 For further reading, the reader is referred to several 
of the numerous excellent reviews that have been published 
recently on the genetics of BMD and hip fracture. 60 " 63 
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Table 3 Loci associated with bone mineral density at the 
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3q25.3l 
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UN7C/DCDC5 
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MARK3 


I4q32.32 


rsl 1623869 




(Continued) 


gQ submit your manuscript | www.dovepress.com 




The Application of Clinical Genetics 2013:6 



Dovepress 



Dovepress 



Clinical impact of genetic discoveries in osteoporosis 



Table 3 (Continued) 



Pathway 


Closest candidate gene 


Locus 


Lead single nucleotide 
polymorphism 


Also associated with 
osteoporotic fracture?* 




MPP7 


lOpI 1.23 


rs3905706 






NTANI 


I6pl3.l 1 


rs4985 1 55 






PKDCC 


2p2l 


rs7584262 






RPS6KA5 


I4q32.l2 


rs 1 286083 


Y 




SLC25AI3 


7q2 1 .3 


rs4727338 


Y 




SMC6 


1 7pl 3.3 


rs479088 1 






SPTBNI 


2pl6.2 


rs4233949 


Y 




STARD3NL 


7pl4.l 


rs69592l2 


Y 




TXNDC3 


7pl4.l 


rs 10226308 






XKR9ILACTB2 


8q 1 3.3 


rs70 1 79 1 4 






ZBTB40 


Ip36.l2 


rs6426749 


Y 



Notes: ^Associated with fracture at P<0.00 1 . Table adapted from Richards et al 55 and Estrada et al. 54 

Abbreviations: Wnt, wingless-related integration site; OPG, osteoprotegerin; RANK, receptor activator of nuclear factor-KB; RANKL, RANK ligand. 



Can genetic discoveries 
help to reduce the burden 
of osteoporotic fracture? 

The last 10 years has been an exciting time for osteopo- 
rosis genetics insofar as GWAS studies carried out during 
this time have led to the discovery of over 60 new loci 



robustly associated with variation in BMD, including some 
subsequently found to be associated with fracture. So how 
might these findings, or future discoveries coming down as 
even larger sample sizes become available, be translatable? 
Will knowledge gleaned from these discoveries improve 
our ability to predict individuals at risk for future fracture, 



Mature osteoclast 



Wnt 



Osteoclast 
precursor 





Bone 
resorption 



Bone 
formation 



Figure I RANK/RANKL/OPG pathway in bone remodeling. The balance between bone formation and resorption is largely regulated by the Wnt pathway (bone formation), 
the RANK {pink symbols)/RANKL (blue symbols) pathway (osteoclast activation), and sclerostin (negative regulation of bone formation). Osteoblasts express the cell surface 
receptors RANKL and Wnt and also secrete a soluble decoy receptor, OPG (green symbols). Wnt protein binds coreceptors Fizzle-Fz and LRP5/6, leading to stabilization 
of p-catenin and its translocation to the nucleus to regulate target genes, resulting in increased bone formation. In the absence of OPG, RANKL on the osteoblast surface 
is available to bind RANK present on osteoclast precursors. Binding of RANK/RANKL leads to osteoclast maturation and resorption of bone. Sclerostin, secreted by 
osteocytes, inhibits Wnt from binding LRP5. 

Abbreviations: RANK, receptor activator of nuclear factor-kappa B; RANKL, receptor activator of nuclear factor-kappa B ligand; OPG, osteoprotegerin; Wnt, wingless- 
related integration site; LRP, low-density lipoprotein receptor protein. 
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allowing initiation of early treatment? Will the discovered 
loci provide novel insights about bone biology and sug- 
gest new therapeutic targets? What new genetics/genomics 
approaches will be applied and what are the prospects for 
their translatability? 

Prediction 

There is good evidence that treating individuals who have 
already experienced a low trauma fracture, or who are 
osteoporotic by virtue of low BMD, with medications to 
either reduce bone resorption or increase bone formation can 
reduce the incidence of subsequent fractures. 64 Unfortunately, 
about one-half of hip fractures occur in women without a 
a pre-existing fracture and who have a BMD that is not in 
the osteoporotic range, 65 and prediction algorithms based 
on BMD and clinical information fractures (eg, DXA and 
FRAX) are imperfect for predicting individuals destined to 
experience a future fracture. An important question therefore 
is whether addition of genetic variants can improve prediction 
accuracy for fracture beyond that which can be determined 
by clinical risk factors alone. 

Based on prospective studies, the discriminatory value of 
prognostic models based on clinical risk factors is moderate 
to good with areas under the receiver operating character- 
istic curve ranging from 0.7 to 0 . 8. 66,67 Tran et al evaluated 
the improvement of these models for fracture prediction that 
would be achieved by adding into the model a variant in the 
collagen I alpha 1 (COL1A1) gene (rs 18000 12) that has been 
consistently associated with fracture risk (reviewed by Tran 
et al). 68 This variant, which marks an Spl binding site in the 
first intron of COLIA1, has a minor allele frequency in Cauca- 
sians of approximately 19%, and individuals homozygous for 
this variant have a 2^1-fold increased risk of fracture. Adding 
this variant to the hip fracture prediction model increased 
the area under the receiver operating curve only modestly 
from 0.86 to 0.88. 68 There was likewise only a very modest 
improvement in reclassification from normal to high risk. In 
the large Genetic Factors of Osteoporosis Consortium meta- 
analysis, a summary genotype score based on genotypes for all 
62 BMD-associated loci had modest ability for prediction of 
osteoporosis (area under the receiver operating characteristic 
curve 0.59), but addition of this score to a model that included 
age and weight increased predictive discrimination only neg- 
ligibly (from area under the receiver operating characteristic 
curve of 0.75 to 0.76). 54 It should not be surprising that the 
marginal increase in discriminatory value achieved by adding 
genotypes at BMD-associated loci is minute given their small 
effect sizes. Moreover, one might expect variants associated 



with BMD to add very little predictive value to models that 
already include BMD as a clinical variable. 

Although single variants have very small effect sizes 
and do not add to prediction, there is a theoretical benefit 
to considering multiple risk variants in aggregate in genetic 
profiling. One problem is that many of the known variants 
associated with fracture risk exert their effects through BMD, 
and thus their additional predictive value would be expected 
to be very marginal. The number of SNPs associated with 
fracture independently of BMD is very small. In simula- 
tion studies, Nguyen and Eisman estimated that 50 SNPs 
with odds ratios of 1.3-1.5 would be needed to provide a 
significant increase in reclassification. 69 However, given the 
complexity of the fracture phenotype, it seems unlikely that 
many associated SNPs will have effects this large. In fact, 
even in the absence of BMD screening, it seems unlikely at 
this point that SNPs will play a meaningful role in fracture 
prediction. 

Novel pathways 

A second area where genetic discoveries could potentially 
be translatable for osteoporosis lies in identifying new 
pathways that could serve as therapeutic targets. Despite 
the growing numbers of SNPs associated with BMD, few 
have led to major new insights about bone biology to date, 
largely because the genes marked by these SNPs have 
either fallen into pathways already known to be related to 
bone biology or the associated SNPs are either intergenic or 
appear to be correlated with genes having functions which 
are not clearly known. To make more sense of GWAS find- 
ings, more comprehensive analyses may be required that 
integrate additional information, such as tissue expression 
of implicated genes and functional annotation of identified 
SNPs and genes. In the large Genetic Factors of Osteoporo- 
sis Consortium meta-analysis, many of the identified SNPs 
mapped to genes in the RANK-RANKL-OPG, mesenchymal 
stem cell differentiation, and Wnt signaling pathways (see 
Estrada et al 54 and Table 3). Integrating expression data to 
GWAS findings can prioritize the most biologically relevant 
candidate genes, as in the GWAS findings carried out by the 
Framingham Study 70 and Genetic Factors of Osteoporosis 
Consortium. 54 In an analysis of 180 SNPs associated with 
height in prior GWAS meta-analyses, Lui et al used rat, 
mouse, and human databases to identify genes expressed in 
growth plates to help localize the causal gene within many 
of the associated loci. 71 

One important application of new discoveries will be in 
the development of medications that more precisely target key 
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mechanistic pathways, such as the RANK-ligand pathway, 
cathepsin K inhibition, and Wnt signaling manipulation. 
Identifying the key players in these pathways, especially the 
ones most amenable to targeting, represents an active area 
of ongoing research. 72 For example, the antiresorptive agent, 
denosumab, is an antibody that inhibits RANKL and has been 
on the market for 3 years, and an sclerostin (SOST) antibody 
has recently completed 5 years of testing and is likely to be 
approved soon. Although these particular targets were identi- 
fied before the advent of G WAS, it is possible that genomic 
approaches will lead to new discoveries in these or other 
pathways that will prove druggable. 

New approaches/future directions 

As for many other complex diseases, the future of genomics 
research in osteoporosis will almost certainly include a thor- 
ough investigation of the role of rare variants discovered by 
deep sequencing of osteoporosis candidate genes as well as 
those discovered through a variety of other sources, includ- 
ing exome sequencing and the 1000 Genomes Project. An 
association study targeted to rare variants will be potentially 
useful for two reasons. First, because GWAS SNP arrays do 
not generally tag variants with minor allele frequencies less 
than 5% well, direct genotyping of these SNPs is the only 
way to assess the low frequency spectrum of allelic variation. 
Second, rare variants located within the exons of genes may 
include some that disrupt gene function, thus potentially 
having large effects on phenotype expression. Even if these 
mutations are rare, their identification can shed important 
insights into the role of that particular gene (and perhaps 
gene pathway) in bone biology. 

Other approaches for identifying genetic variants asso- 
ciated with BMD and/or fracture risk are in their infancy. 
Epigenetic factors that do not modify the DNA sequence 
itself, but rather gene expression, may play important roles. 
The most well known epigenetic mechanisms are DNA 
methylation and histone modifications, which act at the level 
of gene transcription, and miRNAs (microRNAs), which act 
at the post-transcriptional level. Epigenetic marks may drive 
differentiation programs for cell fate, including osteoblastic 
lineage differentiation, and/or may be important regulators of 
bone remodeling during osteoclastogenesis. 73 In fact, recent 
evidence suggests that methylation-dependent mechanisms 
may influence the transcription of RANKL and OPG expres- 
sion. 74 One speculation is that epigenetic marks may be modi- 
fied by nutritional factors, intrauterine growth influences, 
and/or other environmental factors, although at the present 
time it is not clear what these modifying factors are and what 



are the key epigenetic marks most relevant to osteoporosis. 
Resolving these questions has been difficult for a number of 
reasons, including obtaining accurate measurement of the 
relevant exposures, identifying the relevant epigenetic marks 
to measure, and the cell and/or tissue specificity of epigenetic 
marks. These challenges notwithstanding, elucidation of 
epigenetic control of bone regulation may have translational 
implications; for example, some demethylating agents are 
already used to treat neoplastic disorders. 75,76 

There has been very limited research to date into the genetic 
determinants of response to treatment with antiresorptive 
agents used to slow bone loss in osteoporosis. With increased 
knowledge of the molecules and pathways that drive the bal- 
ance between bone formation and bone resorption, one can 
imagine that drug treatments could be refined and optimized 
according to the patient's individual genetic profile to allow 
the most beneficial treatment for each patient. The degree 
to which this can be realized is hard to say, because it may 
depend not only on how individual polymorphisms influence 
response to treatment, but perhaps also on how combinations 
of polymorphisms in aggregate influence response. 

Conclusion 

Over 60 variants identified over the past 7 years have been 
reliably associated with variation in BMD, some of which are 
also associated with fracture risk. Many of these variants map 
to genes in known pathways that are related to bone biology, 
although some map to genes with a relation to bone biology 
that is not readily apparent. A key point is that the effect sizes 
of variants are very small, because they collectively account 
for less than 6% of the variation in BMD. Although predicting 
individuals at risk for future fracture would be of enormous 
value so that treatment can be initiated early, the current panel 
of variants is not useful for prediction of future fracture, given 
that adding them to prediction models that already include 
age, gender, and BMD offers only a negligible improvement 
in predictive value. One hope is that future genetic discoveries 
made through even larger meta-analyses or new discoveries 
based on analyses of rare variants will identify new pathways 
relevant to the pathophysiology of bone metabolism or lead 
to a more granular characterization of known pathways. Such 
knowledge could inspire the development of more targeted 
therapeutics to prevent age-related bone loss. 
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